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Abstract 

We present the exact solution of a model of interacting fermions in any dimen- 
sion with a pure repulsive interaction projecting out a given Cooper channel. 
The solution rests upon the infinite ranged character of the interaction in real 
space, leading to a functional integral that is dominated by a Gaussian term. 
The solution produces strong superconducting enhancements and quasi long 
ranged order in a channel that is not present in the Hamiltonian explicitly, 
but of the form given by arguments from order by projection. 
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There is considerable recent interest |T| in the possibility of models displaying a super- 
conducting behaviour driven solely by repulsive interactions. The search is motivated by 
the High Tc problem, where no obvious,known attractive interaction can account for the 
phenomena; so one believes that repulsive interactions, originating in the Coulomb repulsion 
expressed within the Wannier basis of a few tightbinding bands, are ultimately responsible. 
In one dimension, the usual kind of repulsive interacions generically lead to enhanced SDW 
order, via powerlaw correlations with small decay exponents rather than superconducting 
enhancements. However recent interesting work || shows that under certain conditions, 
superconducting correlations of non trivial symmetry are enhanced. In the physically im- 
portant case of two dimensions the situation is not completely clear in relation to popular 
models, such as the Hubbard or the t-J model. While a treatment within the Random Phase 
Approximation near an Antiferromagnetic instability leads to d-wave superconductivity [[3] , 
one may worry whether the conclusions based on the (weak coupling) approximation are 
valid for strong repulsions. Similar worries exist regarding various versions of the Gauge 
theories || that are in vogue currently. In this context, repulsive models which can be 
solved exactly have an important role to play. A new set of models was introduced recently 
that demonstrate the possibility of enhanced superconducting correlations quite explic- 
itly. These models are generalizations of the Hubbard model, and include a term that is best 
interpreted as a projection operator that excludes a certain Cooper pairing channel from the 
problem. We present here an exact solution of the basic model involving the kinetic energy 
and the projection operator. Our solution is obtained by exploiting a certain feature of the 
interaction within the framework of functional integrals: namely that the projection oper- 
ator is an infinite ranged repulsive operator, and its Hubbard-Stratanovic (HS) linearizing 
field is a single spatially uniform dynamical mode that admits only Gaussian fluctuations in 
the thermodynamic limit .This is in contrast to models with attractive interactions, where 
a new saddle point value of the HS field gets stabilized at low temperatures, and its fluc- 
tuations can be ignored. The repulsive case is much more subtle, and has some similarities 
to the situation found in Fermionic models in infinite dimensions |§ where the relevant HS 
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field is a Grassman variable that factorises the kinetic energy, and has only Gaussian fluctu- 
ations ||. We find that as a consequence the fermi gas develops enhanced correlations in a 
"compromise" pairing channel which is not explicitly present in the Hamiltonian. These are 
precisely of the sort that is expected from arguments from order by projection |5||J . The 
earlier treatments [|S],^| have used a variety of non perturbative techniques, such as rigorous 
inequalities and variational approaches. The exact solution presented here is consistent with 
these, and give us in addition, a clear understanding of the origin of these enhancements 
and of the quasi LRO in terms of the singularities of the two particle scattering amplitude. 
The model is defined by the Hamiltonian 

H = T +U S B ] B . (1) 

Here B = XX(^) ^fc is a Cooper pair operator, b k = C-kiCk] are the pair destruc- 
tion operators, T = ^-k^ka is the kinetic energy corresponding to a band dispersion 
gfc[= — 2 J2a=i cos(k a ) in the nearest neighbour problem ]. may be chosen arbitrar- 

ily. The two cases of interest in 2-dimensions are (i) ((h) = 1 giving rise to extended s-wave 
order, and (ii) ((k) = cos(k x ) — cos(k y ) giving rise to second-neighbour d-wave order. For 
simplicity of presentation we initially focus on £ = 1 and return to the other case later. 

With the above choice of the model, we note that B may be rewritten as J2 r c -l( r ) c t( r )5 
and hence the interaction term may be viewed as an infinite ranged hopping term for doubly 
occupied sites. 

Using the H-S linearization within the standard Grassman variable [ID| formula- 



tion, we write the partition function for this model as the functional integral Z = 
J Dc*Dc Dip* Dip exp— /3^f. The free energy functional ^ is given in terms of the fermi 
fields c, c* and the auxiliary bose fields 0, (p* as 



= [ 0*(r)0(r)rfr-^ f 4 jCT (t)(3 t - ^)c k>tT (r) dr 

-iju s [\f(r)B(r) + £*(t)0(t)) dr, (2) 
v Jo 

where = e k — fi. We use a Fourier series expansion c& )(T (t) = J2 n ex P(^Vi' r )cfc,cr(n), and 



0( r ) = J2n exp(zf2 n r)0 n , where the fermionic frequencies uj n = (2n + l)7r ( 5~ 1 and the bosonic 
frequencies Q n = 2n7if3~ 1 . We denote the Fourier components bk{n) = J2 m C-k,i(rn)ck,-\(n — 
m) and so define B n = J2k bk{ n — !)• In terms of these we can rewrite 



* = E C0n + Efe - «0„)c£ )<y (7l)c fc>ff (n) + l^L/. E(C S n + 4> n B* n ). 

k,n n 

We now trace out the Fermi degrees of freedom and find the reduced free energy func- 
tional ^t/) — + J2(pn(pn — kf,TJ2kTr\og{l — C{k)}. Here \l/o is the noninteracting free 
energy, C is an infinite dimentional matrix defined by its (frequency space) matrix ele- 
ments C n ^ m (k) = U s J2iGo(k,n)Go(k,m)(f)^ n+l+1 (l) n+ i + i in terms of the free Green's function 
G (k,n) = l/(iu n — , and the trace is in the frequency space. 

We now analyse ^ in detail. We can expand 

n 

u 2 . . . . 

Here Uo(ifl n ) = ^ ^fc,™ where itk,n are the polarizations 

7r fc>n = (2/(fc) - l)/(zO n - 2£ fe ) , (3) 

and f(k) is the usual noninteracting Fermi function. The fourth order term is given explicitly 

as 

-p 1 V"^ i^fik) 1)(-Dmi,fc + D m2 ^k)3m 1 +m2,ni+n2 

1 (mi ,ma ,rn,»i2) ~~ "7T n n n n ' 

*■» f, ^- j 'mi,k^m2,k^ni,k^n2,k 

where D m , fc = (iQ m - 2£ fc ) . 

We note that both Ho and T in the above equations are of 0(1) since these are normalized 
sums over momenta. Indeed every term in the expansion has a similar structure and is of 
the same order, namely 0(C). Hence one has the remarkable exact result that the Gaussian 
term dominates the rest of the terms in the thermodynamic limit [JTTJj . Roughly speaking, the 
Gaussian piece gives us the typical size n ~ and so the quadratic ( in U s ) piece is of the 
order 1/C, and likewise the m th term is of the order 1/C^ m ~ x '. Thus in the thermodynamic 



limit, it suffices to keep the Gaussian term and to drop the remaining terms. This leads to 



the following remarkably simple result [12 



From the same arguments, the correlation function of the Cooper pair operators is given by 

< b* k (n)b k i(n) > = 8 k ,k' ^k,n ~ 4 ^k,nVeff{i^n)^k',n (5) 



where the effective interaction 12 



v "' {n) = i + vSo(in.) (6) 

in terms of the (very large) coupling constant V = U S C. 

We next present an alternate derivation of the above results starting from the equations 
of motion, which gives some more insight into them. 

Define (as in ref 0) the set of operators // = J2( e k) l bk, and TJ = XX e fc)'(^-fcj. + n fcT — 1) 
. Clearly I = B and Ti = T. It is easy to ascertain that 

[/,, T] = 2I l+1 , [/,, N] = 2k , [/,, I J] = -T l+m , (7) 

whence, 

[/,, H-fiN} = -2fili + 2I l+1 - U.TJQ. (8) 

Now, we invoke the law of large numbers and argue that in the present problem the operator 
product TiIq can be replaced by C^Iq where C \i[ = < TJ > , the thermodynamic average 
of Ti, which is clearly of order £ [O] . Then the equations of motion for the usual time ordered 



Green functions « If, J TO ' >> reduce to the closed set 

(iQ n + 2fi) « If, I J » = -Cfn +m + 2 « I t+ i;lJ » -CUsM « I ;lJ » (9) 
which can be solved exactly. The solution is given by 

«If,lJ » = -cu l+m + CU t -f— CU m (10) 

1 + U s L,LLo 



where 



,fo (ifi„ + Itf +1 V - 2 &) 

This may be verified using the easily derived recursion relation IL„ +1 = [(zO„ + 2/^)11/ — 

Hi\/2. Furthermore, one of the key features of the models being discussed is that the one 



particle propagators are unrenormalised [14 , so that < n& >= /(&) whence, Hi(iQ n ) = 
(1/jC) J2k ( e k) l/ Xk,n- Then it is easily seen that the results in Eq. (|T0|) are basically the same 
as in Eq. ([5]). Starting from the latter, multiplying by {ek) l {^k') m and summing over k and 
k', we get the former. 

An analysis of the detailed properties of the function V e fj is of crucial importance for the 
rest of our discussion. For a simple model bandstructure, with a constant density of states 
g(e) = 1/2 for —1 < e < 1, we can compute it exactly as 



with 



U (u) = - log( 11 1 { -) + i-9{2 - |2/i + w|)apn(t«;) 



The band extends from —2 + 25 to 2 + 25, where the hole filling 5 = 1 — p = — /i. The 
schematic behaviour of 14// is as follows. At very high frequencies, V e ff ~ V. It has two 
poles at certain large frequencies that essentially dominate its physics . In between these 
poles lies the intermediate frequency range where V e ff is of order unity, which contains the 
branch cut corresponding to the one electron band of states. The location of the poles can be 
found by using a large frequency expansion for the function U (iQ n ), which for |O n | >> W 
( W is the band width) behaves as 

The higher order terms can be verified to be negligible. The poles of the V e ff can then be 
found by solving for the zeroes of the denominator, a quadratic in £ = l/i£l n , given as 
= V" 1 — 8( — 2(|//i| — /i<5)C 2 - The roots are always real corresponding to real frequency 
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poles of V e ff. Hence, (a) at half filling, 5 = 0, the poles are at ±u>o = ±J2V\^i\, and near 
the poles V e f$ ~ ^ 2(0^0) " ^ ne P°^ e fr ec L uenc i es are infinite in the thermodynamic limit, 
but leave behind consequences in the groundstate as we see later, (b) away from but close 
to half filling , 5 << 1, one pole is at U2 = V5 with residue V 2 5, and the other at a large 
negative frequency —u\ = —(2\fii\/5 — 2fi) with residue —uf/5. As 5 — ► 0, the latter poles 
smoothly go over into those of the first case. The contribution from the branch cut does 
not have any particularly simple form, but is not important in the most interesting region 
of the problem, namely 5 ~ 0. 

Using the above properties of V e ff we can compute exactly the instantaneous (expecta- 
tion) values of the the extended s-wave and the s-wave correlation functions : 
a = < A A > with A = [T, B] = -2 £ fc e k b k , and 8 ee I < B*B >. 
Using the main result Eq([|), we find 



a = k B T^[M^n)-V eff (tQ n )Ul(zQ n )]e^ ^ 

n 

fi = ^UV-V eff (M n )}e^ + \ (13) 

The main contributions to the resulting frequency sums are quite easily seen to be tied 
to the pole contributions of V e ff, the branch cut part giving an uninteresting subleading 
contribution. We find at half filling 



Near half filling we find 

"poZe — ^ , JJpo'e — pjJ2£2' ^ °J 

At half filling Eq . (|14|) gives us the quasi long ranged order as well as the correction to ground 
state energy along with their appropriate coefficients, the latter by an integration over the 
coupling constant U s . These answers are in good numerical agreement with exact numerics 
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on a finite sized system Away from half filling we find that the uncertainty principle 
lower bounds are off by a factor of 2, i.e., from Eq.(p3|) a pole = 2aLowerBound- 

The results given above enable us to compute several other response functions exactly. 
For example we find that the charge stiffness or Meissner fraction is non- zero, and in fact 
unchanged from the non interacting value at half filling. We also find at half filling that the 
appropriate order parameter density A(r) has correlations that are novel : < A^(r)A(0) >~ 
-^=. Thus although we do not have LRO of the usual sort, one has a divergent "structure 
function" J df < A f (r)A(0) >. 

We next discuss the important case of £(/c) = cos(k x ) — cos(k y ) in two dimensions. This 
corresponds to suppressing d-wave order at the length scale of nearest neighbours. From the 
uncertainty principle argument of Refs |5]||, it follows that this would lead to enhancement 
again in the d-wave channel, but at a longer length scale, i.e. the resulting A = —2 e kC(k)bk 
should have enhanced correlations. The functional integral solution sketched here bears 
this out exactly. We recover the results in Eq.s([L4],[15]) with the replacements : 5 — > 8 = 
£~ 1 Ek( 2 (k)(l-2f(k)) and fj, -»• /ti = C' 1 £ fe e k ( 2 {k){l -2f(k)). The enhanced correlation 
function of A then diverges at the point where 5 vanishes. So long as the one electron 
dispersion has the bipartite symmetry, one can see that 5 vanishes exactly at half filling. 
However, if the dispersion does not have this symmetry, e.g. by having a second neighbouring 
hopping t', then 5 vanishes at some other density determined by t' , as illustrated in Fig.l. 
The case of t' ~ —.4 is popular in literature, since it leads to a Fermi surface that is 
consistent with that seen in the photo-emission experiments in High T c systems ||15|| , and it 



is an amusing coincidence that the filling 8 ~ .18 is close to the optimum doping . 

Finally, we can show that many of the above features of the infinite range model retain 
their relevance for more realistic models with finite range repulsion. To see this, consider 
(in the s wave case) a more general model of the form H = T + ^ 2 U(p) B'(p)B(p), where 
B{p) = J2k c -fc,i c fc+p,T i s a Cooper pair operator with total momentum p , with U (p) = U s {l c ) d 
for a small set of C/(l c ) d points surrounding the centre of the Brillouin zone. This would 



correspond to repulsive interactions with a long but finite range l c [16|. Then it is not hard 



to see that the results we have discussed above would retain much of their validity (but for 
some differences in details) to leading order in (l// c ), with the replacement U S C — > U s (l c ) d , 
leading to enhanced pairing correlations over a finite range in momentum space. The quasi 
Long Ranged Order would be replaced by an enhancement of the 0(1), like that in the single 
mode model away from half filling. In this case, we see that all Cooper pairs with a finite 
( but small) center of mass momentum are also influenced by the interaction, and thus the 
model is more realistic by way of helping current carrying states. 

In summary, we have found an exact solution for an interesting model of Fermions with 
purely repulsive interactions with infinite range, which may be regarded as a meanfield 
repulsive model. The resulting solution has quasi long ranged order at half filling, as well 
as large unbounded enhancements as one approaches half filling, in the equal time pairing 
correlations. We have also argued that the above methods and results retain their relevance 
even for generalised models where the repulsive interactions have a large but finite range, and 
are therefore more realistic. Finally, it is remarkable that the enhanced pairing correlations 
in these models arise from very high energy poles in the scattering amplitude, not unlike the 
physics of Mott Hubbard systems, where the upper Hubbard band influences the properties 
of carriers in the lower band. 

One of us (BSS) is grateful to the Aspen Center of Physics for support, where a part of 
this paper was written. He thanks P W Anderson and B Shraiman for stimulating discussions 
and constructive suggestions. 
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FIG. 1. The effective filling 5 for d-wave superconductivity against the true filling 5 for 
different values of the second neighbour hopping t'. 
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